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It is shown in previous works that the coupling between two Majorana end states in superconduct- 
ing quantum wires leads to fractional Josephson effect. However, in realistic experimental conditions, 
multiple bands of the wires are occupied and the Majorana end states are accompanied by other 
fermionic end states. This raises the question concerning the robustness of fractional Josephson ef- 
fect in these situations. In this work, we show that the absence of the avoided energy crossing which 
gives rise to the fractional Josephson effect is robust, even when the Majorana fermions are coupled 
with arbitrary strengths to other fermions. Moreover, we calculate the temperature dependence of 
the fractional Josephson current and show that it is suppressed by thermal excitations to the other 
fermion bound states. 



Introduction — Majorana fermions in condensed 
matter systems have drawn much attention re- 
cently because they obey non-Abelian statistics and 
have potential applications in quantum computa- 
tion which are protected from local perturbations 
[112] . A number of systems are proposed to 
host Majorana fermions, including p x + ip y super- 
conductors [3], topological insulator/superconductor 
heterostructures[l], spin-orbit coupled semiconductor- 
superconductor heterostructure [5-9 and metallic thin- 
films in proximity to superconductors jlOj. Majorana 
fermions in these systems manifest themselves through 
several interesting effects. For example, they induce 
cross- Andreev reflections [TT], resonant Andreev reflec- 
tion [12], electron teleportation [13] etc. Moreover, when 
two Majorana end states from two separate supercon- 
ducting quantum wires are brought close to each other 
to form a Josephson junction, they induce fractional 
Josephson effect [T4HT8] . The periodicity of the resulting 
Josephson current is Air in terms of the phase difference 
between the two superconducting quantum wires, instead 
of 27T as in usual Josephson junctions. A schematic pic- 
ture of the Josephson junction is shown in Figflji. 

When two Majorana end states are brought close to 
each other to form the tunnel junction, they form a 
fermion state which may be occupied or empty. These 
two states form the q-bit in proposed quantum computa- 
tion schemes. The fractional Josephson effect is impor- 
tant because its sign depends on the occupation of the 
fermion state and provides a direct way to read out the 
q-bit [HHH]. 

However, in more realistic experimental conditions, 
multiple bands of the quantum wires are occupied and 
the Majorana end states are accompanied by other 
fermion end states. A schematic picture is shown in 
FigjT|3. These fermion end states may couple to the 
Majorana end states strongly. In this case, the effective 
Hamiltonians used in previous works, which describe the 
coupling between two Majorana end states only, do not 
apply. Thus, it is important to understand the robust- 
ness of the fractional Josephson effect when the Majorana 
end states are coupled to other fermions. 

It has been recognized that the presence of other 
fermion states may not affect the function of the Ma- 
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FIG. 1: a) Two ID superconducting wires are brought to- 
gether to form a Josephson junction. The Majorana end 
states at the junction are labeled by vl and rjR respectively. 
These two Majorana end states at the junction couple to 
each other with amplitude T(<f>) to form an Andreev bound 
state. The spectrum of the Andreev bound state is shown 
in Fig[2^. b) When multiple bands are occupied, the Majo- 
rana end states are accompanied by other fermionic states, 
denoted by 7l/_r„ and —jL/Rn (occupied and empty states 
respectively) . 



jorana q-bit, as long as they are localized in the vicinity 
of the Majorana fermions. This is because as long as the 
temperature is much lower than the bulk band gap An, 
thermal excitations of extended states are negligible and 
the total fermion parity, defined as the even or odd oc- 
cupation of the fermions in the vicinity of the Majorana 
fermion, is conserved |19j . However, in the context of the 
fractional Josephson effect, the following questions are 
raised. First, the fractional Josephson effect relies on the 
robustness of the avoided crossing of energy levels formed 
by the Majorana fermions as a function of <fi, the phase 
difference between the superconductors forming the junc- 
tion. Is this avoided crossing robust in the presence of 
strong coupling to many fermion states? Second, what is 
the temperature dependence of the Josephson current? 
Third, while thermal excitation of the fermionic states 
preserves fermion parity, how does the difference of the 
fractional Josephson current which is proposed to mea- 
sure the fermion parity depend on temperature? These 
are the questions we answer in this paper. 

Review of the Single- channel Case — It was first 
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pointed out by Kitaev [TB] that a single channel quan- 
tum wire with a proximity induced p-wave pairing gap 
can have Majorana end states localized at the ends of 
the quantum wire. When two superconducting wires 
with Majorana end states are joined together to form a 
Josephson junction, the two Majorana end states tjl and 
r/R at the junction couple to each other. The effective 
Hamiltonian of the junction can be written as: 

H lDeff = 2iTr lL r m = 2r($^ A - ~), (1) 

where T is the coupling strength and ipA = Vl + ir]R is 
the Andreev bound state formed by the two Majorana 
fermions. Suppose a gauge is chosen such that the su- 
perconducting phase of the wire on the right is fixed, as 
the superconducting phase of the left quantum wire is in- 
creased by 2tt, all the fermions on the left of the junction 
acquire a minus sign such that tjl — > —rjL and t]jt — > r)R. 
Equivalently, the phase dependent can be absorbed into 
r. Thus, F satisfies the condition: 



r(^ + 2rr) = -r(0). 



(2) 



This condition shows that T(4>) is 4-7T periodic and crosses 
zero for certain <fi = <fio. For models in Refs. |6|14)15j . 
T(<fi) = T cos(</>/2), where T is a positive constant. The 
energy spectrum of the junction with T(<f>) — Tq cos(</>/2) 
is shown in Figj2h. 

Moreover, KqTTI shows that T((j>) and —T(4>) are the 
eigenvalues of B.\r> e tt an d V^JO > and |0 > are the cor- 
responding eigenvectors, where ipji\0 >= 0. The energy 
of the occupied states tp\\0 > and the empty state are 
shown by the solid and the dotted lines respectively in 
Fig(2^,. If the fermion parity is conserved at the junc- 
tion, this energy crossing is protected because ^JO > 
and |0 > have different fermion parities and there is no 
transition from one state to the other when <j> equals 2ir. 
Since the energy eigenvalues ±T(<j)) are Air periodic and 
there are no transitions among the states with different 
fermion parity, the Josephson current, which is given by 
I± = ±{2e / h) j^Y [(f)) is also 47r periodic. The observa- 
tion of this 4-7T periodic Josephson current will be strong 
evidence of the existence of Majorana fermions, and the 
sign of the current reveals the fermion parity of the junc- 
tion. 

Multi-channel Case — In more realistic experimental 
situations, multiple transverse sub-bands of the quan- 
tum wires are occupied. It was pointed out in Ref llOl 
that Majorana end states exist in quasi-one dimensional 
wires when an odd number of bands are occupied. These 
Majorana end states are accompanied by fermionic end 
states with nonzero energy and the energy separation of 
the fermionic states is Aq/^/M where A is the bulk 
superconducting gap and M is the number of occupied 
transverse sub-bands. When two quasi-one dimensional 
wires are brought together to form a Josephson junction, 
the Majorana end states are coupled to each other and to 
the nonzero energy fermionic states as well. In this case, 




FIG. 2: a) The energy spectrum of the Andreev bound state 
when the state is occupied (solid line) and empty (dotted 
line). The crossing is protected by fermion parity conser- 
vation, b) The energy spectrum of Hm (solid curves) and 
Pf[/iA/](0)/Pf[/iA/](O) (dashed curve) as a function of <j> are 
plotted. The presence of other fermions induces two more zero 
energy crossings. In this plot, every sign change of Pf[/iAf] is 
accompanied by an energy level crossing at zero energy. For 
the plot, -y L / R , n = 2n, /3 n , n '(0) ~ cos(<?!>/2) and N = 8 are 
chosen in H e ff. 



it is not obvious that the fractional Josephson effect is 
still robust. In this section, we show that the coupling of 
Majorana fermions with other fermionic end states may 
induce an additional even number of zero energy cross- 
ings, i.e. the total number of zero energy crossings is odd 
and the fractional Josephson effect survives, independent 
of the strengths of the couplings. 

The effective Hamiltonian in the multi-band quantum 
wire case can be written as: 



.n=0^ / fl.n=0 + 

E^[27i»(*in**.« -!) + 27fl«(*k, 



R.n 



h.c. 



Here, 



(3) 

'L/R.,n=a = ^\/R n -Q are Majorana end states 
from the left and right of the Josephson junction respec- 
tively. ^L/R,n with ti^O are fermionic states with pos- 
itive energy ^L/R.n- Here, N is the number of fermion 
end states which are within the bulk gap in each quantum 
wire, fi n ,ri (0) are the coupling strengths, n and n' in the 
second sum run from to N . The prime in the second 
sum indicates that the /?oo(0) term is taken out from the 
sum and shown explicitly. As mentioned above, a phase 
change of the fermions on the left by 27T transforms ^ l to 
—^l- This transformation is equivalent to transforming 
fin n' t° —fin,ri- Thus, fi n , n ' satisfies the condition: 



/3 n ,„'0 + 2tt) = -fi n y{^>). 



(4) 



Since the fermionic modes are coupled to the Majo- 
rana modes, it is convenient to use a basis which consists 
of Majorana modes only. Therefore, we write ^l/ro — 

r} L / R o and ^ L /R,n = VL/R,n + iv' L /R,n^ ^L/R.n = VL/R, n - 

ii]' L / R n for the n^O modes, where r\ and rf are Majorana 
fermions. In this Majorana fermion basis, the Hamilto- 
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nian can be written as: 
H M = {TjL,m)hM ( ) = (VL-'in) 



A B 

-B C 



-T 

Tl 

nh 



(5) 

, where rf a = (Va,o,ilaA,Va,i---Va,N,v' a ,N) with a = L or 
R. A and C are (2JV + 1) x (2JV + 1) anti-symmetric Her- 
mitian matrices, which describes the coupling between 
Majorana fermions on the same side of the junction. B 
is a (2iV + l) x (2N + 1) matrix with arbitrary pure imag- 
inary entries, which describe the coupling between Majo- 
rana fermions from opposite sides of the junction. Thus, 
Km is an (47V + 2) x (47V + 2) anti-symmetric Hermitian 
matrix. 

According to Eq|4j the sub-matrices satisfy the condi- 
tions: 

A(6 + 2tt)=A(6), C(6 + 2n) = C(6) 

and B(6 + 2tt) = -B(<j>). {> 

Therefore, changing the phase difference by 2tt is equiv- 
alent to applying a unitary transformation on Km'- 



h M {(t) + 2TT) 



W T h M {<f>)W, 



(7) 



where W = cr z ® I and / is a 27V + 1 dimensional 
identity matrix. Using the fact that Pf[W T tiM{<t>)W] = 
Pf[h M ]Bet[W} and Det[W] = -1, we have: 



Pi[h M (6 + 2n)] = -Pi[h M (</»)], 



(8) 



where Pf denotes the Pfafhan of an antisymmetric ma- 
trix. Since Pf [/im](0) is a continuous function of <fi 
and Pf[/ijv/](0) changes sign when the phase is increased 
by 2n, Pf [Hm] (<t>) must cross zero an odd number of 
times when (f> is advanced by 2tt. From the relation 
Det[/iA/] = Pf[/iM] 2 , vanishing Pf[/ijvf](0) indicates the 
appearance of zero energy modes. 

Moreover, the eigenvalues of Km come in pairs of the 
form {ei(<f>), — £j(</>)} because of particle hole symmetry, 
where i runs from 1 to 27V + 1 and the convention that 
e.j > at 6 = is chosen. As a result: 



2N+1 



(9) 



Suppose zPf[/iM](</>) changes sign at <f> = (f>o, Eq[9] im- 
plies that an odd number of ej changes signs at tfi . Since 
Pi{fiM](4>) crosses zero an odd number of times and each 
crossing is accompanied by an odd number of energy level 
crossings at zero energy, thus, the total number of energy 
level crossings at zero energy must be odd. 

One example which illustrates the relation between en- 
ergy level crossings and the sign of the Pfaffian is shown 
in Fig{2b. It is clear from Figj^ that P{[h M ](2n) = 
— Pf[/ijv7](0) and every sign change of Pf \hM\{4>) is accom- 
panied by a level crossing at zero energy. With an odd 
number of zero energy crossings, a phase change of 2tt 
in a system described by the Hamiltonian Hm does not 



change back to its original state. Instead, a phase change 
of 4-7T is needed. A physical consequence of this 4tt peri- 
odicity is the fractional Josephson effect. Therefore, the 
fractional Josephson effect is robust, even when the Ma- 
jorana fermions are coupled to other fermions strongly. 
This is the first main result of this paper. 

Finite Temperature Current — In order to demonstrate 
the non-Abelian properties of Majorana end states, the 
Majorana end states may undergo braiding operations 
using T-junctions suggested in RefHTl These braiding 
operations may result in changes in the fermion parity 
of the Josephson junction. In this section, we show that 
the Josephson current can be used to detect this change 
in fermion parity at the junction even at finite tempera- 
tures. 

When there is only one single Andreev bound state at 
the junction and when the fermion parity is conserved, 
the Josephson current is 



T + 2e d U\ 



(10) 



where ±eo(</>) are the energy levels of the Andreev bound 
state when the Andreev bound state is occupied or empty 
respectively. In other words, the sign of the Josephson 
current gives the fermion parity of the junction. The 
Josephson current is temperature independent as long 
as thermal excitations to the extended fermionic states 
above the gap are negligible, so that fermion parity is 
conserved. 

In the multi-channel case, if Ao 3> fc^T > Aq/-</M, 
fermion parity is still conserved but the n = state can 
now exchange its occupation with the n ^ fermion 
states. As a result, the Josephson current is temperature 
dependent. In order to calculate the Josephson current, 
we use the relation [20] : 



h 86 ' 



(11) 



where F is the free energy of the Josephson junction. 
Suppose there are TV = \J~M fermionic Andreev bound 
states with energy e„ > e , where eo is the lowest positive 
energy state at 6 = 0, the free energy can be written as: 

1 N 

F e/o = — p lnfQ e^(e^P e/o + e~^P o/e )]. (12) 



i>0 



Here, fermion parity conservation at the Josephson junc- 
tion is imposed such that the free energy with even and 
odd parity are denoted by F e and F a respectively. P e and 
P denote the thermal weights of having an even or odd 
number of e„^o states excited, which can be written as [ 
EI]: 



P. 



1 d 



N 



/o ^ k\ dx k . 

k—even/ odd i—1 



ri( 



(13) 



The finite temperature current can be obtained from 
Eq |ll| with two contributions to the fractional Josephson 
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current. The Josephson current originated from the eo 
state: 

^ w= T{ e^; o+e -^p o ; e } # 60(ft (14) 

and an extra contribution originated from the e„^o 
states: 

^'yl 'lt'-'f 1 ' (15) 

To extract the fractional Josephson current from a 
background 2tt periodic Josephson current, we define the 
quantity: 

h/o{4>) = l[lTe/o(4>) - lTe/o((t> + 2tt)] 

= l[I e /o{<t>) ~ Ie/ott + 27T)] + |[^ /o (0 " ^ /o (0 + 2tt)] 

(16) 

where 7r e and /to denote the total Josephson current 
when the fermion parity of the junction is even and odd 
respectively. There are two contributions to I e / {<P)- The 
contributions from the eo state and the e n ^o states are 
given by I Fe/o = |[7 e/o (0) - I e / (<j> + 2tt)] and I' Fe/o = 

|[7g/ c (<^) — 7g/ o (0+27r)], respectively. A nonzero value of 

I(4>) signals the presence of fractional Josephson effect. 
Moreover, using the fact that eo(4>) = —eo(cf) + 2ir) and 
e n #o(<?f>) = e n #o(0 + 2tt), we have: 

U = -L (17) 

This is the second main result of our paper, which shows 
that the fractional Josephson current continues to depend 
on the fermion parity of the junction when the Majorana 
fermions are coupled to other fermionic modes and at 
finite temperatures. 

In the multi-band superconducting wire case, it is rea- 
sonable to consider the situation in which the two Ma- 
jorana fermions ry^o an d Vro are weakly coupled to each 
other such that eo *C 5e, where <5e = e n +i — e„ is the 
energy level spacing of the fermionic states. 5e is in 
the order of Ao/y/M when the coupling between the 
fermionic end states are weak [10] - In the regime where 
5e > k B T « e , we have P e sa 1 and P Q sa e ~ 2tlf} < 1. In 
this case, Ifc/o ~ g^ £ o- On the other hand, I' Fe i ~ 
± 2£(- e -2 £o /3 _ e 2c /3) e -2 £l ^^ e ^ Therefore, the I' F term 

is exponentially suppressed. The fractional Josephson 
current in this regime is given by the zero temperature 



result of Eq[l0] multiplied by a temperature dependent 
prefactor which is exponentially small when k B T > e\. 
The fractional Josephson current and the maximal am- 
plitude of the current is plotted in Fi gj3] A simple mod el 

with e (</>) = r o cos(0) and e„ = \J (nA Q /^/M) 2 + e^j is 
assumed. Fig(3] shows that the maximal amplitude of 
the fractional current decreases exponentially as temper- 
ature increases. However, at temperature k B T ss <5e, the 
fractional Josephson current is of order Tq. Assuming 



Se = 2r = 0.02meV, we have I - InA at T w 0.2K. 




<P k B T/r 



FIG. 3: a) I e and I verse <f> at Se = 2Fo and at temperatures 
ksT = lFo and fcsT = 3Fq. Solid and dotted curves denote 
I e and I respectively. I e / a are in units of j[-Fo- b) The 
maximal amplitude of the fractional Josephson current I c at 
Se — 2Fq and Se — 3To respectively. 



Interestingly, in the regime where ej > > 5e, the 
Ipe/o terms are suppressed and the I Fe / terms domi- 
nate. As a result, the (^-dependent fractional Josephson 
current is very different from the ones in Figfja. This 
regime can be reached when the quasi-one dimensional 
quantum wires become very wide. The details of this 
work will be reported elsewhere. 

In conclusion, we show that the absence of the avoided 
energy crossing, which gives rise to fractional Joseph- 
son effect, is robust even when Majorana fermions at the 
Josephson junction are coupled with arbitrary strengths 
to other fermion states. The fractional Josephson cur- 
rent can be used to measure the fermion parity of the 
junction at finite temperatures, as long as the tempera- 
ture does not greatly exceed the excitation energy to the 
other fermions. 
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